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A doublet of cosmic-ray events with primary energies > 10 20 eV 
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The Telescope Array Collaboration has observed a cosmic-ray event with estimated primary energy of 
1.38 x 10 20 eV whose arrival direction coincides [9], given the angular resolution of 1.5°, with that of an event 
with estimated primary energy of 1.23 x 10 20 eV observed by the Pierre Auger Observatory. The total number 
of events with energies > 10 20 eV published by both experiments is six. I estimate the statistical significance 
of the doublet, which is rather weak, and point out that the arrival directions of events in the doublet coincide 
with the Galactic X-ray source Aql X-l. 



Despite decades of intense studies, including those 
by recent huge experiments, sources of ultra-high-energy 
cosmic rays (UHECRs) remain unknown. For the pri- 
mary cosmic-ray particles with energies of order 10 20 eV 
or higher, quite simple astrophysical arguments re- 
strict the number of potential accelerating astrophysical 
sources drastically [TJ [5] . At the same time, the Greizen 
[3], Zatsepin and Kuzmin [4] (GZK) effect shortens the 
mean free path of protons and nuclei with those high 
energies considerably, putting an additional constraint 
that sources of these events should be nearby. This 
logic has been supported by the recent observations of 
the flux suppression at high energies consistent with 
the GZK predictions by the High-Resolution Fly's Eye 
(HiRes) [S], Pierre Auger Observatory (PAO) [5] and 
Telescope Array (TA) [7] experiments. The observa- 
tion of the suppression does not mean however that 
no "super-GZK" events are observed. Two largest and 
|most modern arrays of surface detectors have published 
coordinates of three events each [5J H] with energies 
E > 10 20 eV. It is these six events which will be pri- 
marily concerned in this note. 

For the "sub-GZK" (E~6x 10 19 eV) events, early 
PAO data suggested a weak correlation of cosmic-ray ar- 
rival directions with positions of nearby active galactic 
nuclei (AGN) [TU] which might be interpreted as an indi- 
cation to acceleration of UHECRs in these astrophysical 
objects. This interpretation has been criticised on the 
basis of numerous arguments, see e.g. [TTJ [T2J H3] ; it has 
not been supported by the data of HiRes [2] and TA [5] 
though it has been supported by the Yakutsk data [15] . 
A subsequent publication of the Pierre Auger collabo- 
ration 8 , based on enlarged statistics, demonstrated a 
much weaker effect. However, the highest-energy events 
with E > 10 20 eV did not correlate with AGN even m 
the data set with the strongest signal. 



Exp. 


Date 


E, EeV 


RA 


DEC 


PAO 


09.10.2008 


123 


287.7° 


+1.4° 


TA 


28.02.2011 


138 


288.5° 


-0.0° 



Table [TJ Details of the two events with coinciding ar- 
rival directions: the experiment name; date; energy in 
units of EeV= 10 18 eV; equatorial coordinates. 
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Figure [TJ The sky map with arrival directions of three 
PAO events with E > 10 20 eV (diamonds) and three 
TA events with E > 10 20 eV (boxes). The Hammer 
projection, equatorial coordinates. 



As it has been pointed out in recent Ref. [5], where 
coordinates of TA events have been made public for 
the first time, the arrival direction of one TA event 
with E > 10 20 eV coincides, within the experimental 
precision, with that of a PAO event of the similar en- 
ergy. Details of the two events are given in Table [TJ 
for convenience. The sky map with all six events with 
E > 10 20 eV is presented in Fig. [TJ The appearence of 
the doublet is psychologically surprising because the two 
experiments are located in different hemispheres and see 
different parts of the sky with a moderate overlap in the 
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Figure [2] The sky map with arrival directions of the 
two events in the doublet: the PAO event (diamond) 
and the TA event (box). With a 68% probability, the 
true arrival directions are inside the corresponding cir- 
cles. The star denotes the position of Aql X-l; no other 
strong X-ray or gamma-ray sources are seen nearby. 

equatorial region. The zoom of the skymap, with error 
circles of events, is presented in Fig. [5] 

To estimate the statistical significance of this dou- 
blet we follow the usual procedure described in Ref. [16] 
(see also [UJ1IB]). We assume the isotropic distribution 
of arrival directions, account for direction-dependent ex- 
perimental exposure and simulate a sufficient number 
of Monte- Carlo sets of events, then count how often the 
same or larger number of doublets happens as a fluctu- 
ation of the isotropic distribution. For this purpose, a 
doublet is denned as a pair of arrival directions sepa- 
rated by not more that v^o where the angular resolu- 
tion 6» « 1.5° for both PAO and TA. 

For three PAO and three TA events (E > 10 20 eV), 
the P-value calculated in this way is P 3.7 x 1CP 3 . 
This value may be interpreted as an estimate of the 
probability to have one doublet anywhere in the com- 
bined data set as a fluctuation of the isotropic distribu- 
tion of arrival directions. However, this interpretation 
should be taken with care because the choice of the en- 
ergy threshold, 10 20 eV, is somewhat arbitrary. 

We see that the statistical significance of the doublet 
is not that impressive. Many three-sigma effects have 
come and gone in cosmic-ray physics. Nevertheless, it 
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Figure [3] The X-ray light curve of Aql X-l (quick-look 
results provided by the ASM/RXTE team [24]). Verti- 
cal lines with arrows denote the arrival times of cosmic 



rays. 



is tempting to speculate about the potential source of 
the particles. The region of interest is located close 
to the Galactic plane, in the zone of avoidance where 
not many extragalactic objects are identified. How- 
ever, strong active galaxies which might accelerate par- 
ticles up to ultra-high energies are expected to be X-ray 
and/or gamma-ray sources visible through the dust ob- 
scuration at this location (Galactic coordinates of the 
center of the doublet are I « 35.9°, b w -4.3°). In 
representative catalogs of active galactic nuclei (Veron- 
Cetti and Veron [12]), gamma-ray sources (2FGL [2D] ) 
and X-ray sources (ROSAT bright source catalog [2T]). 
there is only one bright source within a few degrees of 
this location, a low-mass X-ray binary star Aql X-l, 
just in the middle of the doublet (see Fig. [5]). There are 
no active galaxies nor other identified ROSAT X-ray or 
FERMI-LAT gamma-ray sources around. 

Aql X-l, the brightest X-ray source in the Aquilla 
constellation, is an X-ray millisecond pulsar in a binary 
system (see e.g. Ref. [22] for discussion and references). 
The system is located at the distance of kpc from 

the Earth [23] . It experiences quasi-periodic outbursts 
each 300 days roughly (see the X-ray light curve in 
Fig. [3]) • Though the object is one of only twelve known 
Galactic accretion-powered millisecond pulsars |25] and 
is well studied, it does not appear very exotic. It is sin- 
gled out of this dosen only by a relatively large mass of 
the companion in the binary system, M > 0.45M©, and 
the correspondingly large orbital period of ~ 19 h. The 
estimated magnetic field on the neutron-star surface is 
- (1 ... 5) x 10 8 G [2S] [37]. On the basis of X-ray tim- 
ing and spectral properties, this object is classified as 
"atoll" (see e.g. Ref. [28] for a more detailed discussion 
of classifications). Accretion in these sources may have 
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similarities with accretion on stellar-mass black holes 
[29]. 

An extragalactic E ~ 10 20 eV proton arriving from 
the direction we consider should be deflected by the 
magnetic field of the Milky Way by — (2... 4)°, de- 
pending on the held model which is not known precisely. 
However, a hypothetical particle coming from Aql X-l 
would be deflected by a much smaller amount because 
the source is considerably closer to the Earth than the 
Galcatic Center is. Assuming charge one and the mean 
Galactic magnetic field in the disk of ~ 1 /iG, one ob- 
tains a rough estimate for the deflection of ~ 1.1°. This 
deflection would correspond, for a proton, to the time 
delay of ~ yr, thus making it not surprising that the 
arrival times of the events do not coincide with par- 
ticularly interesting moments in the life of the would-be 
source, cf. Fig. [3] (unless neutral primaries are assumed). 
We note that currently, neither PAO nor TA is able to 
determine the primary particle type of a particular air 
shower detected by the surface array. 

In general, a wide belief that cosmic rays with E > 
10 19 eV are of extragalactic origin is based on a few rea- 
sonable arguments. First, these energetic particles are 
not confined by the Galactic magnetic field and, assum- 
ing similar fields exist in other galaxies, are not confined 
anywhere. Second, the arrival directions of these cosmic 
rays are (almost) isotropic on large angular scales, while 
the distribution of any kind of Galactic sources on the 
sky is anisotropic. Third, there is a lack of Galactic 
objects where sufficient conditions for acceleration of 
particles to those energies are satisfied. Nevertheless, 
some proposals for Galactic sources are being discussed 
(see e.g. Refs. [30l EH [32] ) . 

The first two arguments might be overcome if in the 
Milky Way there are only a few sources of cosmic rays. 
Then, part of the observed events come from these few 
Galactic sources and the rest comes from similar sources 
in nearby galaxies (thus explaining weak hints of correla- 
tion with the local distribution of matter at the highest 
energies). However, Galactic sources should then dom- 
inate the flux at high energies (cf. e.g. Ref. [33 ) thus 
producing either the Galactic anisotropy or a concen- 
tration of arrival directions towards particular sources. 
A price to pay for not seeing this in data is the fine 
tuning which is however not excluded. 

As for the third argument, acceleration of UHECRs 
in pulsars has been proposed long ago in Ref. [31] (see 
Ref. [35] for a different recent proposal). In general, it 
is difficult to overcome radiative energy losses in pul- 
sar magnetosphere [21 135); however, in the regime where 
the losses are dominated by the curvature radiation and 
the electric and magneic fields are parallel in a very 



long tube ("linear accelerator"), the energy- loss limits 
may be relaxed [2] [35]. Another problem with pul- 
sars is the screening of the accelerating potential gap 
when electron-positron pairs are created; however, one 
expects that some pulsars are "pair-starved" [36] [37] 
and may accelerate particles to higher energies. It is 
presently unclear whether these conditions are satisfied 
in Aql X-l. 

To summarize, we observe an interesting coincidence 
of the arrival directions of two out of six cosmic-ray 
particles with estimated energies in excess of 10 20 eV 
observed by the modern surface-detector arrays, the 
Pierre Auger Observatory and the Telescope Array ex- 
periment. The probability that one or more doublets 
occur in the data set by a fluctuation of the isotropic 
flux is about 3 x 10~ 3 . The error circle of the arrival 
directions of coinciding cosmic rays includes a bright 
accretion-powered millisecond pulsar in the X-ray bi- 
nary Aql X-l. Future studies are required to support or 
reject the conjecture that Galactic sources are able to 
produce super-GZK cosmic rays and whether this par- 
ticular source possesses physical conditions allowing for 
particle acceleration to that high energy. 
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sian Federation NS-5590.2012.2 and by the "Dynasty" 
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